Rotational properties of dipolar Bose-Einstein condensates confined in anisotropic 

harmonic potentials 
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We study the rotational properties of a dipolar Bose-Einstein condensate confined in a quasi-two- 
dimensional anisotropic trap, for an arbitrary orientation of the dipoles with respect to their plane 
of motion. Within the mean-field approximation we find that the lowest-energy state of the system 
depends strongly on the relative strength between the dipolar and the contact interactions, as well 
as on the size and the orientation of the dipoles, and the size and the orientation of the deformation 
of the trapping potential. 
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I. INTRODUCTION 

The rotational properties of Bose-Einstein condensates 
have been studied extensively over the last decades, be- 
ginning with the well-known experiments performed on 
hquid Helium. Bose-Einstein-condensed gases of ultra- 
cold atoms have provided us with an ideal system for 
the study of this problem In addition to the fact 

that they are dilute, many of the system parameters are 
tunable externally, making it possible to explore experi- 
mentally many novel properties. So far most of the ex- 
periments have been performed in isotropic as 
well as anisotropic [12] harmonic potentials. More re- 
cently, other trapping potentials such as, for example, 
anharmonic [l^ and toroidal [ijj ones, have also been 
addressed. 

While the initial experiments considered atoms with 
interactions that could be modelled via an isotropic and 
local, i.e., contact, potential, more recently the behav- 
ior of atoms with an electric or magnetic dipole mo- 
ment has been investigated jisl , [Tg} . The physics of these 
gases is rather different because the dipolar interaction 
is anisotropic, non-local, and finally it can be both at- 
tractive - when the dipoles are placed head-to-tail - and 
repulsive - when they are placed side-to-side. Initially 
dipolar effects were investigated experimentally in ^^Cr 
atoms [Tt'I. More recently dipolar effects have also been 
investigated experimentally in condensates of '^^K and 
^Li [ill atoms, in spinor condensates [l^, in ultracold 
gases of Dy atoms ^20] , and in ultracold polar molecules 
|2ll [22j . These studies have focused on different aspects 
of the dipolar interaction such as, e.g., its anisotropy Ts*], 
or on its attractive character and the consequent possible 
collapse of the gas 22] . 

The rotational properties of Bose-Einstein condensates 
are described e.g. in the review articles [InSj, as well as 
in Refs. (23 - l30| . In particular, anisotropic trapping po- 
tentials have been considered in Refs. |26l - l3C| . Dipolar 
gases have also been extensively addressed theoretically, 
see for example Refs. [isl lia [3ll - [39| . and studies on dipo- 
lar Bose-Einstein condensates under rotation have shown 



that the vortex properties of these systems are strongly 
affected by the anisotropic character of the dipole-dipole 
interaction and by the relative strength between the dipo- 
lar and the contact interactions j34l. IsgI, IstI. |40|. |41| . The 
relative orientation of the dipoles with respect to the con- 
fining potential also affects the critical rotation frequency 
for the nucleation of a vortex state in the condensate, low- 
ering (raising) its value with respect to the non-dipolar 
case when the trap is elongated in the direction perpen- 
dicular (parallel) to the orientation of the dipoles [35|| . 

Some of the above studies have employed a fully three- 
dimensional approach to solve the Gross-Pitaevskii equa- 
tion (35I Isrj. whereas other works have examined quasi- 
one- or quasi- two-dimensional systems [34, 36, 38, 39). In 
most of these papers, the orientation of the dipoles has 
been limited to only two possibilities: perpendicular, or 
parallel to their plane of motion. Very recently, the more 
general and interesting case of an arbitrary orientation 
of the dipole moment of the atoms has been investigated 
for non-rotating gases [H, [s^ [42| . 

Here, we study the rotational properties of a quasi-two- 
dimensional dipolar Bose-Einstein condensate confined in 
a rotating anisotropic harmonic trapping potential, with 
an arbitrary orientation of the dipoles. Since the lat- 
ter can be controlled by means of external electric or 
magnetic fields, this is also very relevant from the exper- 
imental point of view. As compared to the case where 
the gas is trapped in an axially-symmetric trap and inter- 
acts with a contact potential, here the symmetry is bro- 
ken in two separate and independent ways: Through the 
anisotropy of the trap, and also through the directional 
dependence of the dipolar interaction. The combination 
of these two symmetry-breaking mechanisms gives rise to 
interesting vortex configurations, which are investigated 
in this article. 

In what follows we first present our model in Sec. II. In 
Sec. HI we investigate the rotational properties of the sys- 
tem as a function of the trap anisotropy, the strength of 
the dipolar interaction, and the orientation of the dipole 
moment of the atoms with respect to their plane of mo- 
tion. Finally, in Sec. IV we summarize our results and 
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FIG. 1: (Color online) Schematic illustration of the quasi- 
two-dimensional elliptical confining potential. A very tight 
potential is assumed along the z axis (not shown in the plot). 
The atoms have a dipole moment in the xz plane, at some 
direction that forms an angle 6 with the x axis, due to the 
action of an external electric (or magnetic) field E. 



conclude. 



II. MODEL AND METHODOLOGY 

Let us consider a Bose-Einstein condensate of atoms 
of mass M confined by some external potential, which is 
harmonic in all directions, 



(1) 



Here is the frequency of the potential along the z axis, 
which we assume to be very tight, while in the xy plane 
the potential is taken to be anisotropic, 

^ \M{u:lx'+^y). (2) 

In the above expression rj^ — {x,y), and uj^. = 
idoV^ + 4A, LOy = cjoVl — 4^ are the trap frequencies 
along the x and y axes, respectively. Here A is a dimcn- 
sionless parameter determining both the direction and 
the strength of the deformation of the trap. 

The assumption of a very tight potential along the z 
axis, huz being much larger than any other energy scale 
in the problem, allows us to treat the problem as quasi- 
two-dimensional. With this assumption, the degrees of 
freedom along the z axis are frozen and the system occu- 
pies only the lowest-energy eigenstate (t>o{z) of the poten- 
tial V{z) = Muj'^z'^/2. As a result, the order parameter 
^'(r) separates to \I'(r) = ^{r±)(j)o{z). We also assume 
that the atoms have a non-zero (electric or magnetic) 
dipole moment, and that all the dipoles are aligned to 
some external (electric or magnetic) field, which is in the 
xz plane, and forms an angle Q with the x axis, as illus- 
trated schematically in Fig.[T] In three dimensions, the 
dipolar interaction between two atoms separated by the 



vector r is given by [16j 

Vddir) = D- 



1 - 3 cos^ 



(3) 



In the above expression, — (P / (47reo) when the atoms 
have an electric dipole moment d, where cq is the permit- 
tivity of the vacuum. When the atoms have a magnetic 
moment /x, = hqh"^ / (At:), where /io is the permeabil- 
ity of the vacuum. Also, Ord is the angle between the 
polarization direction and the relative position vector r 
of the two dipoles. 

As shown in Refs. [s^, [s^, in order to derive an ef- 
fective potential for the two-dimensional problem, one 
assumes the Gaussian profile of the harmonic oscillator 
in the z direction 4>a{z)- Then, integrating Vdd{^) over 
the z coordinates, 



\M^)\''Vdd{v-T')\M^')\'dzdz'. (4) 



Working with cylindrical polar coordinates r — (rj^ , 
14ff(r±) = -J=^|(2 + 2w)Kq(w/2) ~ 2wKi{w/2) 



cos^ e 



-(3 + 2w)Ka{w/2) + (1 + 2w)Ki{w/2) 



-2 cos^ 9 cos^ ( 



-wKq{w/2) + {w~- l)Ki{w/2) 



.(5) 



Here Uz = 
rection, w 



y^Ti/mujz is the oscillator length in the z di- 
= r\/2al, and Koiw) and Ki{w) are the 
zero-order and first-order modified Bessel functions of 
the second kind. It can be checked easily that for the 
particular cases = 0° and 90^ Eq. ([5]) reduces to the 
expressions found in Refs. [s^, HBl- The assumed quasi- 
two-dimensional behavior of the system allows us to de- 
rive the above expression for the dipolar interaction, and 
also allows us to consider any arbitrary orientation of the 
dipoles. 

In addition to the dipolar interaction, we also consider 
the usual contact potential between the atoms, given by 
Mnt(r) = UoS{r), where Uq = Anh^a/M. Here a is the 
scattering length for zero-energy elastic atom-atom colli- 
sions. Again, integrating over the profile along the z axis, 
the corresponding effective two-dimensional contact po- 
tential becomes Vs{r±) = gS(r±), where g is given by 
g = UoJ \Mz)\^dz = Uo/{V2^az). 

Assuming also that the trap rotates around the z axis 
with some angular frequency fi, the Gross-Pitaevskii 
equation for the order parameter ip{r±) takes the fol- 
lowing, nonlocal form 



2M 



+ l^(r±) + Fdip(ri) + ff|V^(ri)P - QLz 



= A'V'(r-L), (6) 

where Lz is the operator of the angular momentum along 
the z axis, /i is the chemical potential, and 



^dip(ri) = I Kfr(ri - r^) |^(rl)|' dr'^ 



(7) 



3 



is the dipolar interaction potential. We have treated this 
term making use of the convolution theorem and fast- 
Fourier-transform techniques, combined with the intro- 
duction of a cutoff at small distances, where Vcg{r±) di- 
verges, as done in Rcf. 32]. 

In order to solve Eq. ^ we have employed a fourth- 
order split-step Fourier method within an imaginary-time 
propagation approach [i^. Thus, starting with a reason- 
able initial state we propagate it in imaginary time until 
we reach a steady state, after a sufficiently large number 
of time steps. 



edd=0 edd=0.31 6^^=0.38 edd=0.54 £^,,=0.67 




III. RESULTS - ROTATIONAL PROPERTIES 
OF A DRIVEN GAS 

Since there is a large number of parameters in the 
present problem, we fix ujz and ujq throughout this study, 
choosing uJz/ujo = 100. 

One useful dimensionless quantity is the ratio between 
the interaction energy per particle due to the contact po- 
tential Eg and the oscillator energy TiutQ. For a cloud 
of homogeneous dens ity of r adius equal to the oscil- 
lator length oq = ^/h/MujQ the interaction energy is 
Eg = gN/ (TTag), where N is the atom number. The ratio 
Ea/hujQ is thus equal to gNM/ (nh^), and in what follows 
we choose it equal to IOO/tt. 

Another dimensionless quantity that is convenient to 
introduce is the ratio between the "dipolar length" add, 
and the s— wave scattering length a, which we denote as 
Edd = CLdd/a Defining add = MD^/{3h'^), then Sdd = 



4:TTD^/3Ua- Since Uq ~ \/2TTgaz, Edd is also equal to 
{^/8^^ /3)D^ / {gaz) . This quantity gives roughly the ratio 
between the expectation value of the dipolar energy £^dip 
and Eg. Therefore, with the above choice of parameters, 
and for the specific values of D that we choose below, we 
have the following hierarchy of energy scales: 



hojz > Es> Edip > hujo 



(8) 



It should be noted that the stability of the system is al- 
ready a non-trivial question [H, , as the dipole-dipole 
interaction is partly attractive and one has to be cau- 
tious, since the system may be unstable against collapse. 
We have observed this collapse in our calculations, al- 
though not in the range of parameters corresponding to 
the figures shown below. 



A. Vortex configurations as a function of the 
dipole strength and the dipole orientation 

Let us first study the rotational properties of the sys- 
tem as a function of the orientation of the dipoles of the 
atoms with respect to their plane of motion, and of the 
dipolar strength. In Fig. 2 we show the two-dimensional 
atom density for = 0, 15, 30, 45, and 90 degrees, and 
Edd = 0, 0.31, 0.38, 0.54, and 0.67. We also choose a 



FIG. 2: (Color online) Two-dimensional atom density as a 
function of the strength of the dipolar interaction, for Edd = 
0, 0.31, 0.38, 0.54, and 0.67, and of the orientation of the 
dipoles, for = 0, 15, 30, 60, and 90 degrees, with the non- 
dipolar case (sdd = 0) also shown as a reference. The trap 
deformation is A = —0.03, and the angular velocity of the 
trap is Q/iuo = 3/4. The inset illustrates the coordinate axes 
and the external field. 



moderate value oi A = —0.03, which corresponds to a 
trapping potential that is more tight along the y axis and 
finally ft/uJo is set equal to 3/4. From Fig. 2 one can see 
that the way the vortex structure changes as a function of 
Edd depends sensitively on the orientation of the dipoles. 
For example, when the polarization is along the x axis 
(O = 0°) the attractive part of the dipolar interaction 
gives rise to a well-known self-induced squeezing of the 
density in the y direction [s^] that reduces progressively 
the number of vortices as Edd increases. 

Thus, the four vortices present in the non-dipolar case 
(edd = 0) and for Edd = 0.31, get reduced to three and to 
one for Edd = 0.38 and 0.54, respectively, with a struc- 
ture that clearly reflects the symmetry imposed by the 
external fleld along the x axis. Finally, when the dipo- 
lar strength exceeds a certain value (Edd — 0.67), the 
width of the gas becomes too small to accommodate any 
vortex state. A similar effect has been found in non- 
dipolar condensates which are conflned in anisotropic po- 
tentials, where there is a critical deformation of the trap 
bey ond which all the vortices are expelled from the cloud 

When the polarization angle becomes Q = 15° we find 
a similar, but less pronounced disappearance of the vor- 
tex lattice, as in this case the attractive part of the dipo- 
lar interaction becomes less important. For Q — 30°, 
however, the repulsive and the attractive parts of the 
dipole interaction become comparable, and as a result 
the vortex structure is not affected by the increase in Edd 
within the range that we have considered. 

When the angle exceeds 60° the behavior of the system 
changes. In this case, the repulsive part of the interaction 
becomes dominant, giving rise to a stronger net repulsion. 
As a consequence, the number of vortices increases with 
increasing Edd- This effect is maximal when the dipoles 
are polarized along the z axis {Q = 90°), in which case 
the dipolar interaction is purely repulsive and isotropic. 
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A=-0.1 A=-0.05 A=0 A=0.05 A=0.1 




FIG. 3: (Color online) Two-dimensional atom density as a 
function of the ellipticity of the trap for A — —0.1, —0.05, 0, 
0.05 and 0.1, and of the orientation of the dipoles, for = 0, 
15, 30, 60 and 90 degrees. The angular velocity of the trap 
is fl/ujo = 3/4, and Edd is fixed to 0.67. The coordinate axes 
correspond to the same as in Fig. 2. 



B. Vortex configurations as a function of the 
ellipticity of the trap and the dipole orientation 

We have also studied the rotational behavior of the 
system as a function of the trap deformation A and of 
the angle Q, for a fixed dipolar strength Edd = 0.67 and 
a fixed rotational frequency of the trap, fi/ojo = 3/4, as 
shown in Fig. 3. In particular, we considered the cases 
of strong and moderate deformation of the trap in the 
y axis {A — —0.1 and A = —0.05, respectively), the 
non-deformed case {A = 0), and of strong and moderate 
deformation along the x direction {A = 0.1 and A — 0.05, 
respectively). When the dipoles are polarized in-plane, 
along the x axis {Q = 0°), the squeezing of the density in 
the y direction reported in the previous subsection due 
to the attractive part of the dipolar interaction prevents 
the formation of vortices, not only for A < 0, but also 
in the non-deformed case {A — 0). However, when the 
trap is squeezed in the direction of the polarization of 
the dipoles, the confinement competes with the dipolar 
interaction, forcing the cloud to expand along the y axis, 
and thus allowing the formation of vortices in the sys- 
tem, with their number increasing as the trap gets more 
deformed. 

The same qualitative behavior is found when the polar- 
ization angle 6 is slightly increased. Indeed, for Q = 15° 
the main difference with the previous case is that the 
cloud has a vortex state for A = and A = —0.05. Sim- 
ilarly, when the angle becomes 30 degrees, the tendency 
of reducing (increasing) the number of vortices when the 
trap is squeezed in the y (x) axis is still clear, although in 
this case the vortex lattice does not disappear completely, 
even for the largest deformation that we have considered. 

This trend is however no longer observed for 9 = 60°. 
In this case, the dipolar interaction is mostly repulsive 
and the number of vortices increases with respect to the 
non-deformed case, regardless of the sign of the defor- 
mation. When the latter is moderate, the vortex struc- 
ture is the same for A = ±0.05, and differences are only 



present in the case of strong deformation. Finally, when 
the dipoles are polarized along the vertical direction, i.e., 
when Q — 90°, the dipolar interaction is isotropic and 
purely repulsive and as a consequence the densities are 
axially symmetric for A = and have a mirror symmetry 
with respect the interchange of the x and y axes in the 
deformed cases. 



C. Angular momentum versus the rotational 
frequency of the trap 

Finally, we have studied the expectation value of the 
angular momentum of the system in the z direction, (L^), 
as a function of the rotation frequency il, for different 
values of Q and A, and for a fixed strength of the dipolar 
interaction Sdd = 0.67, as shown in Fig. 4. 

The different cases shown in Fig. 4 share some com- 
mon characteristics. First of all, they are consistent 
with the divergence of the angular momentum when 
57 minjoj^:, ujy}, i.e. when D, LUy — woVl — 4: A since 
A is chosen to be positive. Also, as ft/uiQ approaches the 
limiting value y/l — AA, the curves become more smooth, 
as the number of vortices increases and the system be- 
comes more classical, resembling in this limit solid-body 
rotation. We also observe in Fig. 4 that for some given 
Q, and A, as 8 decreases, the angular momentum of the 
gas decreases as well. 

In addition, the results in Fig. 4 are consistent with 
the intuitive expectation that the more the system is dis- 
torted from axial symmetry (i.e., as A increases, or as Q 
decreases) , the more it resists in accommodating a vortex 
state , and as a result the critical frequency Q.c neces- 
sary for vortex nucleation increases. This effect is most 
clearly seen in Fig. 5, which shows the critical value of VI 
for the nucleation of the first vortex state as a function 
of <d. This increase of ilc is associated with a distortion 
of the axial symmetry of the cloud, which may carry an- 
gular momentum even when it is vortex-free, as required 
by the irrotational nature of the velocity field. 

In most cases shown in Fig. 4 there is a discontinuous 
increase of (L^) with increasing f2, which is associated 
with the progressive entry of vortex states into the gas, 
very much as in the case of non-dipolar atoms confined 
in an isotropic trapping potential [25|. One interesting 
exception is shown in the top panel of Fig. 4, where for 
= 0° (i.e., when the dipolar interaction breaks the ax- 
ial symmetry in the most extreme way) the angular mo- 
mentum increases continuously from zero. In this range 
of il/ojQ the density of the cloud is elongated along the 
direction of the polarizing field, breaking the axial sym- 
metry of the potential, and it is also vortex-free. Beyond 
a critical rotational frequency {ri/ujQ « 0.82), however, 
the angular momentum shows a sudden increase, which 
is associated with three vortices entering the cloud. On 
the contrary, for the other values oi A = 0.03 and 0.1 
that we have considered, there is a discontinuous transi- 
tion to a state with one vortex, as seen in the central and 
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FIG. 4: (Color online) Expectation value of the angular mo- 
mentum of the gas as a function of the rotational frequency 
of the trap, for y4 = (top), A = 0.03 (middle) and A = 0.1 
(bottom), for Sdd = 0.67, and = 0, 30, 60, and 90 degrees. 
The dashed vertical lines in the lower two curves correspond 
to the value Q/luq — \/l — 4A. 

bottom panels of Fig. 4, where (L^) jumps from zero to 
unity. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have investigated the stationary states 
of a dipolar Bose-Einstein condensate that is confined 
in a rotating, elliptical, quasi-two-dimensional trapping 
potential. In particular, we examined this problem as a 
function of the orientation of the dipole moment of the 
atoms, of the relative strength between the dipolar and 
the contact interactions, and of the deformation of the 
trap. 

While a repulsive contact potential between the atoms 
tends to spread the gas uniformly, the trap deformation 
and the anisotropic character of the dipolar interaction 
favor the breaking of the axial symmetry. These two 
mechanisms give rise to interesting configurations, which 
we have investigated here. For example, at a fixed trap 
deformation, we observe either the disappearance of the 
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" 15 30 45 60 75 90 
(degrees) 

FIG. 5: (Color online) Critical rotational frequency of the 
trap for the nucleation of the first vortex as a function of Q, 
for A = -0.03, and for Sdd = 0.04, 0.17, and 0.67. 



vortex lattice, or nucleation of vortices as the dipolar 
strength increases, depending on the orientation of the 
dipoles. A similar behavior is found if we fix the strength 
of the dipole interaction and we vary the orientation of 
the dipole moments, or the deformation of the trap. 

Another interesting observation is that even in the 
case of an axially-symmetric trapping potential, rela- 
tively weak values of the dipolar interaction energy (as 
compared to the energy due to the contact potential), 
suffice to give rise to rather significant effects due to the 
breaking of the axial symmetry. An analogous result has 
been found in the case of dipole-free atoms that rotate 
in an asymmetric trapping potential, where even weak 
deviations from axial symmetry have a substantial effect 

One rather general and intuitive observation that re- 
sults from the present study is that the more the gas is 
distorted from axial symmetry (via either the deforma- 
tion of the trap, or via the tilt of the angle of the dipoles 
from the direction perpendicularly to their plane of mo- 
tion), the more it expels the vortices, with the critical 
frequency for vortex nucleation being increased j24j . 

The strength of the dipolar interaction affects the crit- 
ical frequency for the formation of vortices in the cloud 
in a non-monotonic way. When the dipoles are oriented 
along their plane of motion, or close to this direction, 
the rule mentioned in the previous paragraph applies 
and, and as a result, the critical rotational frequency for 
vortex nucleation increases with the dipolar interaction 
strength. On the other hand, when the dipoles are either 
perpendicular to their plane of motion, or close to that 
direction, the critical rotational frequency for vortex nu- 
cleation decreases as the dipolar interaction increases. As 
noticed in Ref. [sl] this is due to the fact that when, for 
example, O = 90°, the dipolar interaction is purely repul- 
sive and isotropic, and as a result the average atom den- 
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sity becomes lower. While the dipolar interaction affects 
the large-scale properties of the cloud (i.e., its radius), it 
also leaves its small-scale properties (i.e., the coherence 
length) unaffected, which is reflected in the decrease of 

Our study gives just a flavor of the richness of the phys- 
ical effects that result mainly from the anisotropic nature 
of the dipolar interaction, but also from the deformation 
of the trapping potential. In particular, the sensitivity 
of the behavior of the gas on the orientation of the ex- 
ternal polarizing field (and thus of the dipole moment of 
the atoms), which may be tuned easily, is remarkable. 



The experimental investigation of these effects appears 
therefore worthwhile. 
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